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spaced. 



PACS numbers: 04.50.Kd, 04.70.-s, 04.25.Nx 

I. INTRODUCTION 

Some decades ago Regge and Wheeler began a pio- 
neering study of a small perturbation in the background 
of a black hole in order to get information of the stabil- 
ity of this object [I], a problem that was continued by 
Zerilli 2\ . The oscillations found in these studies are not 
normal modes due to the emission of gravitational waves; 
thus, their frequencies are complex and, as a result, the 
oscillations are damped. 

The terminology quasinormal mode (QNM) and quasi- 
normal frequency (QNF), aiming to name these new 
modes and their frequencies, was first pointed out by 
Vishveshwara [3] and Press [I]. Although initially stud- 
ied in black hole backgrounds, the concept of QNM ap- 
plies to a much broader class of systems. The QNMs 
of black holes were first numerically calculated by Chan- 
drasekhar and Detweiler [5] showing that the amplitude 
is dominated by a ringing characteristic signal at interme- 
diate times. The QNMs are independent of the particular 
initial perturbation that excited them. The frequencies 
and damping times of the oscillations depend only on the 
parameters of the black hole and are, therefore, the "foot- 
prints" of this structure. Soon, the connection of QNMs 
to astrophysics was established by noting that their ex- 
istence can lead to the detection of black holes through 
the observation of the gravitational wave spectrum. The 
interest in QNMs has motivated the development of nu- 
merical and analytical techniques for their computation 
(see [5HSj for a review). Also, the study of the quasinor- 
mal spectrum gives information about the stability as- 
pects of black hole solutions using probe classical matter 
fields (scalar, electromagnetic, spinorial) evolving in the 
geometry without backreacting on the spacetime back- 
ground. Much has been done in that direction, not only 
in four dimensions [5] [10] , but also in two [11] , and in 
more than four [T2] . 

Aside from the study of the stability of the solutions, 



the QNFs are important in the context of the gauge- 
gravity correspondence, whose most celebrated example 
is the duality between the type IIA-B string theory in 
AdSs x S 5 spacetime and the four-dimensional supersym- 
metric Yang-Mills theory [13] . Such a correspondence can 
be generalized for those cases in which there is an event 
horizon in the gravity side. In this case the Hawking tem- 
perature of the black hole is related to the temperature of 
a thermal field theory defined at the anti-de Sitter (AdS) 
boundary. Also, as a consequence of the correspondence, 
the quasinormal spectrum corresponds to the poles of 
thermal Green functions |14j . more precisely, the inverse 
of the imaginary part of the fundamental quasinormal 
frequency can be interpreted as the dual field theory re- 
laxation time [15] . 

Another interesting application of QNMs appears in 
the context of black hole thermodynamics. Some decades 
ago Bekenstein [IB] suggested that the horizon area of a 
black hole must be quantized, so that the area spectrum 
has the form A n — ^ynh, with 7 a dimensionless con- 
stant to be determined. The first proposal to calculate 
this constant through QNMs was made by Hod [T7] . Ac- 
cordingly, the real part of the asymptotic quasinormal 
mode can be seen as a transition frequency in the semi- 
classical limit, and its quantum emission causes a change 
in the mass of the black hole, which is related to the 
area. In this way, the constant 7 for a Schwarzschild 
black hole was determined as 7 = 4 In 3. Later, Kun- 
statter |18j repeated the calculation quantizing the adia- 
batic invariant I = J dE/uj(E) via the Bohr-Sommerfeld 
quantization and using the real part of the QNF as the 
vibrational frequency. The result was exactly the same as 
Hod's. However, recently Maggiore [19] pointed out that 
QNMs should be described as damped harmonic oscilla- 
tors, thus, the imaginary part of the QNF should not be 
neglected, and the proper physical frequency is the mod- 
ule of the entire QNF. Moreover, when considering the 
quantization of the adiabatic invariant, the frequency to 
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be considered is that corresponding to a transition be- 
tween two neighboring quantum levels. With this iden- 
tification, the constant 7 for a Schwarzschild black hole 
becomes 7 = 87r, a result that coincides with the value 
calculated by other methods [5U]. The consequences of 
Hod's and Maggiore's proposals were promptly studied 
in several spacetimes [3TJ [55] . 

In this paper, we are interested in the study of the 
stability of the z — 3, three-dimensional Lifshitz black 
hole found in the context of the so-called new massive 
gravity (NMG) [53]. Moreover, as an application of our 
QNM results we aim to calculate the area spectrum of 
this black hole. 

NMG is a novel parity-preserving, unitary |24j , power- 
counting super-renormalizable |25j . three-dimensional 
model describing the propagation of a massive positive- 
energy graviton with two polarization states of helicity 
±2 in a Minkowski vacuum, whose linearized version is 
equivalent to the Pauli-Fierz theory for a massive spin-2 
field in three dimensions. The action of NMG consists of 
a "wrong sign" Einstein-Hilbert term plus a quadratic 
curvature term given by a precise combination of the 
Ricci tensor and the curvature scalar, which introduces a 
mass parameter [26] . As with other massive gravity the- 
ories, NMG also admits black hole-type solutions with 
several asymptotics and additional parameters [371 HI] • 
Even though this last feature could challenge the usual 
Einstein-Hilbert gravity, it is seen that the definition of 
mass in this new type of black holes is a conserved charge 
computed from a combination of the black hole parame- 
ters, which satisfies the first law of thermodynamics. A 
study of QNMs in these static new type of black holes 
has been performed in [29j . 

The black holes we take into account for our study are 
asymptotically Lifshitz, i.e., they exhibit the anisotropic 
scale invariance, t — > X z t, x — > Ax, where z is the dy- 
namical critical exponent. Specifically, we deal with the 
solutions found for the particular case of z = 3 and a pre- 
cise value of the mass parameter [23 . The general class 
of these solutions are important in the context of gauge- 
gravity duality 30, 31 and were also investigated in other 
background theories [3"5H3"5] . No stability study of black 
holes with Lifshitz asymptotics in three dimensions in the 
scenario of NMG has been performed yet. We aim to give 
some contribution to this issue by considering the QNF 
of scalar and spinorial matter fields in the probe limit, 
i.e, there are no backreaction effects upon the asymp- 
totically Lifshitz black hole metric. Spinor fields have 
been extensively studied in general relativity 36] [37] , and 
their quasinormal frequencies have also been considered 

[55] EH] 00]- 

The paper is organized as follows. In section 2 we in- 
troduce the Lifshitz black holes and discuss their causal 
structure. Sections 3 and 4 are dedicated to the study 
of stability under scalar and spinorial perturbations with 
special emphasis on the massless spinor for the latter. In 
section 5 we present the numerical analysis for both kinds 
of perturbations showing the QNMs and the correspond- 



ing QNFs computed in each case. Section 6 is devoted to 
the calculation of the area spectrum of these black holes 
as an application of our quasinormal spectrum. Finally, 
we discuss our results and conclude in section 7. 



II. LIFSHITZ BLACK HOLES IN THREE 
DIMENSIONS 

In this section we review the black hole solutions we 
will consider within this paper, and we comment some of 
their features. 

The NMG theory [55] is defined by the (2 + 1)- 
dimcnsional action, 
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where m is the so-called "relative" mass parameter, and A 
is the three-dimensional cosmological constant. Defining 
the dimensionless parameters, y — m 2 l 2 and w = XI 2 , it 
is found that the theory exhibits special properties at the 
points y = ±1/2. When looking for black hole solutions 
with Lifshitz asymptotics, it is precisely at the point y = 
— 1/2, w = —13/2, with Lifshitz scaling z — 3, where the 
field equations turn out to be solved by [23] 



a 1 J„2 , .2,-2 



ds z = -a(r) — dr + —dr z + r 
H A 



where 
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with M an integration constant and I 2 = — Also, 
the NMG admits as a solution, the well-known Banados- 
Teitelboim-Zanelli (BTZ) black hole with the dynamical 
critical exponent z — 1. As we shall see below in more de- 
tail, this metric ([2| exhibits a regular single event horizon 
located at r = r + = lyM and a spacetime singularity at 
r = 0. Besides, the surface r = r + acts as a one-way 
membrane for physical objects as we can see from the 
norm of a vector x normal to a given surface s. Since s 
has to be null in order to be a one-way membrane, the 
norm of x must be null as well, i.e., g rr = 0, which occurs 
at r = r+. 

From the behavior of the Kretschmann invariant for 
the metric ([2]), 



R^XaR^ Xa = ~J4^ [«4 - 48j +r 2 + 91r 4 ] 



we see that, for 
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and for r — > 

R^XoR^ -> oo . (7) 

Thus, the black hole solution has a genuine spacetime 
singularity at the origin r — and an event horizon at 
r = r+. Nevertheless, to see if the singularity is timelike, 
spacelike, or null we have to construct the Penrose-Carter 
diagram. First of all, we must remove the coordinate 
singularity at r = r + . Rewriting the metric in terms of 
null coordinates (U, V) 

U = e r ^ t+r '\ V = -e-^-^ . (8) 

where is the tortoise coordinate shown in the next 
section, we get 



the behavior of the Kretschmann invariant. Moreover, 
it is light-like and covered by a regular event horizon at 
r = r + . 

III. SCALAR PERTURBATION 

In this section, we analyze the behavior of a scalar field 
perturbation in the background of a three-dimensional 
Lifshitz black hole. 

The scalar field obeys the Klein-Gordon equation, 

□$ = -La A f (V^99 MN d N ) $ = m 2 $ , (11) 

where m is the mass of the field $. Performing the de- 
composition 



ds 2 = 



4 \r + 



'dUdV . 



(9) 



which is manifestly regular at r = r+. 

Finally, to construct the Penrose-Carter diagram 
(Figfl]) we use the following set of null coordinates 



T = U + V, X = U-V, 
with U — arctan(f) and V — arctan(V). 




(10) 



FIG. 1. Penrose-Carter diagram for the Lifshitz black hole. 
The singularity at r = is light-like and covered by a regular 
event horizon r+. 

From this diagram we see that the spacetime singu- 
larity is located at r = 0, as previously observed from 



*(t,r,# = *(t,r) e** 
The Klein-Gordon equation takes the form, 
Ml 2 ' 
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Ml 2 



* = 0.(13) 



Even though this equation has an analytical solution, 
as we will see in what follows, it is also useful to check the 
numerical results. With this goal we further decompose 
\& = X(t,r*)/y/f, where the tortoise coordinate is 
given by 



= l 4 



1 ( r 

— — — - arccoth — —= 



1 

Ml 2 r 



(14) 



In this way the Klein-Gordon equation reduces to 



-d 2 X + d 2 ^X = V(r)X , (15) 
where V(r) is the scalar effective potential given by 
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Now let us come back to the issue of finding an exact 
solution for Eq.(13). We set the time dependence of the 



field fy(t, r*) as R[r)e lult and redefine the radial coordi 



nate as r = l\M/y. Thus, Eq.(13) turns out to be 
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R = 0, 



(17) 



whose solution is given in terms of Hcun confluent functions, 

R(y) = C iy 2+a (l - y 2 f /2 HeunC (0, a, 



/3 2 a 2 



4 ' 4 
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+C 2 y 2 - a (l~y 2 f /2 HeunC [Q.-a.jS 



' 4 ' 4 
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(18) 
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where Ci and Ci are integration constants, while a = Ci = 0. In order to apply the boundary condition of 
\/4 + m 2 ? 2 , and /3 = —ilu/M 3 / 2 . in-going waves at the horizon we use the following con- 

Imposing the Dirichlet condition at infinity we set nection formula |41j . 



HeunC(0, 6, c, d, e, z) 



Cl r(l - b)T(c) 
T(l + c + k)T(-b-k) 

C2 r(i - 6)r(-c) 



r(i 



s)T(-b- 



HeunC(0, c, 6, — d, e + d, 1 — z) + 

-(1 - z)" c HeunC(0, -c, 6, -d, e + d, 1 - z) 



(19) 



This formula connects a solution around the singular reg- 
ular point z = to the corresponding solution about the 
singular regular point z = 1 of the confluent Heun equa- 
tion given by 



and e is related to e as 



be 
~2~ 



Thus, near y = 1 Eq.(18) can be written as 



(23) 



z(z-l)i/" + [(6+l)(z-l) + (c+l)z]H' + ((iz-e) J ff = 0. 

(20) 

The parameters k and s are obtained from 



fc 2 + (fc + c+l)fc-e + d/2 = 0, (21) 
s 2 + (6-c+l)s-e + rf/2 = 0, (22) 



tf(z/^i)«£i(i-y 2 ) 



0/2 r(i + a)r(/3) 
r(a - fe)r(i + /3 + fc) 



+^-y^ r(a-.)r(i -/* + .)' (24) 

with £j as constants. As we are looking for quasinor- 
mal frequencies with negative imaginary parts, which 
give stable solutions, we find that for (3 < we need 
r(l + /3 + fc) — > oo. Thus, the quasinormal frequencies 
are 



u> = 2i- 



M 3 / 2 



I Q 2 

+ 2N + ^A + m 2 l 2 - J 7 +-m 2 l 2 + — + (3 + 6N) ^4 + m 2 l 2 + 6iV(iV + 1) 



(25) 



where iV is a positive integer. The imaginary part of the While the asymptotic frequency (N — » oo) is given by 
fundamental frequency (N = 0) is negative provided that 

r- M 3 ' 2 

luoo = -2(V6 - 2)i— — N < 0. (27) 

Thus, since the imaginary part of the quasinormal fre- 
3 k 2 r — r — quencies is negative provided that the parameters respect 



/ 3 ^2 - 

y 7 + -m 2 l 2 + — + 3V4 + m 2 l 2 > 1 + \J A + m 2 l 2 



the relation ( 26 1 , we can conclude that the model is sta- 



(26) ble under scalar perturbations. 
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IV. SPINORIAL PERTURBATION 

In this section, we are going to consider a spinorial field 
as a perturbation in the spacetime given by the three- 
dimensional Lifshitz black hole. We analyze the covariant 
Dirac equation for a two component spinor field VP with 
mass fi s . This equation is given by 
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(«) 



e (a) 'V^-M s * = 0, 



(28) 



where Greek indices refer to spacetime coordinates 
it,r,<f>), and the Latin indices enclosed in parentheses 
describe the flat tangent space in which the triad basis 
G/ a \ is defined. The spinor covariant derivative V p is 
given by 



V, = ^ + k (a)W [7(a), 7(6)] 



where w, 



(a)(b) 



is the spin connection, which can be writ- 



ten in terms of the triad e, as 

(a) 



(30) 



where ^ are the metric connections. The 7^ denotes 
the usual flat gamma matrices, which can be taken in 
terms of the Pauli ones. In this work we take 7^ = £02, 
7' 1 ' = o~ 1 , and 7^) = 03. 

We can write the triad basis for the metric (2 1 as 
follows: 



(a) _ y/ajrjA (a) (a) _ r , (a) 



e 2 (a) =r<5 2 (a) , 



(31) 



(29) and the metric connections, 



1 01 



dr 



In 



a(r)A 



1/2' 



r 1 
1 n 



dr 



In 
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A 
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2r^dr 



a(r)A 



A 

r 



r 2 

1 12 



(32) 



r 



Using these quantities it is straightforward to write At this point we are able to write the Dirac equation 
down the expressions for spin connection components. for the two component spinor 
In the present case, we have only two nonvanishing com- 



ponents, 



(o)(i) _ 1 d ( a(r)A \ (1)(2) 



(34) 



2yfa{r)dr \ r 



33) which turns to be the set of coupled differential equations 



* r ^ -r *v A „ T i i 
d t t> 2 + 5 r * 2 + -d^i + - 



y^)A 
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r r 



a(r)'A A' 
+ 



a(r)r r\J~K_ 

a(r)'A A' 

— V 1 

a(r)r r %/A 



*2-M,*i =0, 
*i - A« s *2 = 0. 



(35) 
(36) 



r 



In order to simplify our problem, we redefine 4"i and Thus, we can rewrite Eqs.([35[) and (36) as 
vp 2 as 



#1 = i [a(r)A] 1/4 e - 4Wt+im ^$ + (r), 
* 2 = [a W A] 1/4 e- it * +im *$_(r) , 



(37) 



and the tortoise coordinate as in the scalar case ( 14 ) 



d 
dr* 



Ayfajrj d 



dr 



(38) 



J air) A , . , , 

t V V (39) 



\/a(r)A , . , 

d Tt $ + + = i v v 2 ; (m + ifi s r) , (40) 

where to = zto. 

Furthermore, we define a new function 0, a new rescal- 
ing for the spinorial components R± , and a new tortoise 
coordinate r* through the expressions, 

6* = arctan(^l^) , 

TO 
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FIG. 2. Decay of scalar field with mass m = 1 and I = 1 for 
different values of black hole mass M. 



<t>± = e 



±i0/2 



R±(r) : 



— arctan(^*— ) . 
2ui rh 



In this way Eqs.(|39[) and (40 1 become 



± iu) R ± =WR T , 
where W is the so-called superpotential, 

iy/a(r)A (to 2 + /i 5 . 2 r 2 ) 3 / 2 



W 



9 / - n 9 9\ /J. s m-\/a(r) A 

r A \m A 4- W*%T ) H 2 



(41) 



(42) 



Notice that when o(r) = 1, Eq.([42j) reduces to the BTZ 
superpotential [42 . 



Finally, letting X± = i?+ ± i?_ we have 

(0?. + w 2 )A± = U±A±, 
where V± are the superpartner potentials, 

dW 



V±=W Z ± 



dr* 



(43) 



(44) 



which in the case of a massless spinor (jj, s = 0) reduces 
to 



, m 2 M mM i-z 
V±= ( — T —^r 2 



Ml' 



to 2to i—z 



Ml- 



(45) 



NUMERICAL ANALYSIS 



In this section, we numerically solve EqsJl5| and (43 1, 



which correspond to the scalar and massless spinorial 
perturbations, respectively. Although in the scalar case 
we found an analytical solution and the corresponding 



r 



QNF, our motivation to perform the numerical analysis 
is to verify the applicability of certain numerical meth- 
ods in asymptotically Lifshitz spacetimes. In particular, 
it would be interesting to check if the Horowitz-Hubeny 
method [T3] works well when finding the QNF. 



Using the finite difference method, we define ip(r*,t) — tp(— j'Ar*, lAt) = "4>j,i, V(r(r*)) = V{— iAr*) = Vj, and 
rewrite Eqs.(15) and (43) as 



tpj,l+i - 2i/)jj + ipj+i.i - tyjj + il>j-i,l 



At 2 



Ar 2 



which can be rearranged as 



At 



1>j,i+i = -tjj-i + ^ {i>j+U + $j-hl) + ( 2 - 2— - At% ) ^ . 



+ 0(At 2 ) + O(Arl) = , 



At 2 



Ar 2 



(46) 



(47) 



The initial conditions ip(r*,0) = fo{r*) and ^(r*,0) = solution is stable if 
fi(r^) define the values of ^pj t i for I — and I = 1 and ^ 2 ^2 
we use Eq. (47) to obtain the values of ijjj t i for I > 1. 1 Vmax < 1 > 



Ar 2 



(48) 



At j = we impose Dirichlct boundary conditions since 
V(r. ) tends to infinity as r, tends to zero. The numerical where y^ = ^ is the largegt value of y. in Qm domail 
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This condition is verified in all cases studied here. 



•* 2 , we obtain 



Now we are going to analyze the potential for the scalar 
case. By rewriting Eq.(16) in terms of a new variable 



V(r) = 



2 /2 

ml 



2 /2 

m I 



K 2 l 2 

Zh 







K 2 l 2 \ 2 ] 


j ZhZ + 


e- 


Zh ) ""I 









(49) 



where = ri. The parable in brackets tends to infinity 
as long as (| + m 2 / 2 ) > 0, which is consistent with the 
Breitenlohner-Freedman-type bound for the present case. 
The roots of this polynomial potential are given by 



zo 


= 0, 




= z h 


Z- 


= z h 



K l 
Zfc 



m 2 l 2 



(50) 
(51) 




FIG. 3. Decay of scalar field with mass m = 1 and I = 1 
varying the azimuthal parameter k . 



If m 2 l 2 > — 1, we see that z_ < z + . Thus, going back 
to the original variable r, the roots of the potential are 
r = with double multiplicity, r = ^J~zZ and r = 77, 
(r = — ^fzZ and r = — are excluded as r > 0). Then, 
since is the biggest root and V(r) tends to 00 when 
r tends to 00, the potential is positive-definite in the 
region (r^,oo). Therefore, the quasinormal modes for 
m 2 l 2 > — 1 are necessarily stable [T5] . 

The numerical results regarding the decay of the scalar 
field are shown in Figs[2][3j and the comparison between 
the numerical and analytical results is displayed in Figj4] 
Our results reinforce the conclusion already found ana- 
lytically; the z — 3 Lifshitz black hole is stable under 




FIG. 4. Imaginary part of scalar quasinormal frequencies. We 
display the results obtained using different approaches. 



scalar perturbations. Moreover, according to FiggJ the 
numerical results have a very good agreement with the 
analytical calculation. 

Figure [4] shows that the Horowitz-Hubeny method 
gives unreliable results. In 48J, it is argued that the 
frequencies do not converge as required by the method, 
and that may be explained by ill-conditioned polynomi- 
als. However, in this work, the frequencies converge, but 
they do not agree with the analytic expression and with 
the results from finite difference method. In [33], the 
authors find cases where this method does not work ei- 
ther, and they do so by comparing the results with other 
methods. For instance, they find that the method is un- 
reliable for dimensions bigger than 6. Even in the original 
work [15] the method is unreliable for small black holes, 
and there is no clear explanation for this limitation. In 
our case the asymptotic behavior of the black hole under 
study might play an important role in the convergence 
of the method. Nevertheless, a general criteria for the 
convergence of the Horowitz-Hubeny method remains an 
open question. 



In the case of the massless spine-rial perturbation the and their derivative turns to be 



superpartner potentials ( 45 1 can be written as 
1 



V± = 



(ml) 2 r 2 (r 2 - r\) ± (ml)r 2 (2r 2 - r\)J: 



r 2 — r3_ 



(52) 



Vl = 4 < (ml) 2 r(2r 2 - r\) ± (ml) 



2r(r 2 - r\)Jr 2 -r\ + 



2r 2 — rj 



(53) 



We can see that outside the event horizon V+ is 
positive-definite if ml > 0, and liirv^oo V + (r) = — oo 
if ml < 0. Whereas V- is positive-definite if ml < 0, and 
lim^oo V-(r) = —oo if ml > 0. Moreover, we notice 
that if ml — 0, we have a free-particle case. The decay- 
ing behavior of the massless spinor is given in Figs|5][6] 
Thus, we conclude that the z = 3 Lifshitz black hole is 
stable under massless spinorial perturbations. 




FIG. 5. Decay of massless spinor with I = 1 and black hole 
mass M = 1.0 for different values of the azimuthal parameter 
m. 



VI. AREA SPECTRUM 

One of the applications of our results for the quasinor- 
mal frequencies is the relation they have with the area 
spectrum of a black hole. According to Maggiore [19], 
the proper physical frequency of the damped harmonic 
oscillator equivalent to the black hole quasinormal mode 



ir 4 





m=1 . 




m=3 — 




m=7 






m=10 - 











FIG. 6. Decay of massless spinor with I = 1 and black hole 
mass M = 1.5 for different values of the azimuthal parameter 
rn. 



is given by 



(54) 



where ojr and u)j stand for the real and imaginary part 
of the asymptotic QNF, respectively. Thus, using (27) 
we have 



lj p = 2(V6 — 2) — - — N . 



(55) 



According to Myung et al. [43], the Arnowitt-Deser- 
Misner (ADM) mass of the Lifshitz black hole we are 
studying is given by 



M = 



M 2 



(56) 



Applying Maggiore's method, we calculate the adiabatic 



9 



invariant / as 



I = 



dM 
Aw 



M 
Auj 



dM . 



(57) 



where Aw is the change of proper frequency between two 
neighboring modes, i.e., 



Alu = 2(V6 - 2 



Thus, 



I = 



IM 1 ' 2 
(V6-2) 



(58) 



(59) 



which is quantized via Bohr-Sommcrfcld quantization in 
the semiclassical limit. Recalling that the horizon area 
of the black hole is given by A — 2irr + , with r + = l^/~M, 
and using (59), we arrive at the result, 



A = 27r(v / 6- 2)nh. 



(60) 



with n an integer number. Therefore, we see that the 
horizon area for the z = 3 Lifshitz black hole is quantized 
and equally spaced. 

This result would not be expected for a theory contain- 
ing higher order curvature corrections since, in general, 
black hole solutions in such theories do not have a propor- 
tional relation between their entropy and area, and con- 
sequently, both of them (if any) should not be quantized 
with an equally spaced spectrum for large quantum num- 
bers [HI 2H US] • However, it was already demonstrated 
that the z — 3, three-dimensional Lifshitz black hole has 
an entropy proportional to its horizon area |43U46j . Thus, 
our result (60) also states that the entropy should be 



quantized with a spectrum evenly spaced. Nevertheless, 
we should stress that a generalization of this result for 
Lifshitz black holes should wait for the calculation of the 
area spectra of other black holes of such a type. Solely 
these studies can give a definite answer on this subject. 



VII. CONCLUDING REMARKS 

We have studied the stability of the three-dimensional 
Lifshitz black hole under scalar and spinorial pertur- 
bations in the probe limit through the computation of 
quasinormal modes. In addition, we have found the event 
horizon area quantization as an application of the results 
for quasinormal modes using Maggiore's prescription. 

Regarding the stability, we have not found unstable 
quasinormal modes in the range of parameters that we 
have considered; all the frequencies have a negative imag- 
inary part indicating that the modes are damped and 
thus, the perturbations decay, leaving the system stable 
against this particular sort of probe fields . 

In the case of a scalar probe field, such results totally 
agree with the analytical expressions for the quasinormal 
frequencies; they show a very large imaginary part and 



a very small real part. These modes are almost purely 
imaginary. We have implemented two different numeri- 
cal methods in order to obtain the quasinormal frequen- 
cies and modes: the finite difference and the Horowitz- 
Hubeny methods. The former allows us to obtain the 
temporal behavior of the fields showing all the stages of 
the decay, while the latter gives only the frequencies val- 
ues. As explained in section [Vj the Horowitz-Hubeny 
method failed in the calculation of the scalar frequencies 
as it can be observed in Fig(3] On the contrary, the fi- 
nite difference method has a very good agreement with 
the analytical expression ( |25[ ). Apart from the numeri- 
cal factor, the asymptotic scalar frequency found in the 
present work is the same as the one calculated in the 
hydrodynamic limit of the scalar perturbations in the 
context of gauge-gravity duality [4"T]. 

Regarding the spinorial perturbation, our numerical 
results show that the probe massless spinor decays and, 
thus, the 2 = 3 Lifshitz black hole is stable also under 
spinorial perturbations. 

As a by-product we also obtained the area spectrum 
of this black hole by means of the application of Mag- 
giore's method using our results for the scalar asymptotic 
quasinormal frequencies. Equation (60) shows that the 
horizon area is quantized and equally spaced. Further- 
more, in light of the conclusions shown in [?3J ESL the 
corresponding entropy should also have an evenly spaced 
spectrum. 

Finally, although we have demonstrated the stability 
of the z = 3, three-dimensional Lifshitz black hole under 
scalar and spinor perturbations, we should stress that the 
definite answer on stability should come from the gravita- 
tional perturbations, in particular, from the tensor part 
of the metric perturbation. It is well known that Einstein 
gravity in three dimensions has no propagating degrees 
of freedom, however, the massive versions of the theory, 
e.g., NMG, allow the propagation of gravitational waves. 
Albeit this subject deserves further study, the calculation 
of metric perturbations is a formidable task that is out 
of the scope of the present paper. The analysis is not 
dead easy because the perturbation equation is a fourth 
order differential equation. Thus, some other techniques 
need to be used together with the traditional QNM anal- 
ysis (50|. This study will be discussed elsewhere. 
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